Fitting CMB data with cosmic strings and inflation 
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We perform a multiparameter likelihood analysis to compare measurements of the cosmic mi- 
crowave background (CMB) power spectra with predictions from models involving cosmic strings. 
Adding strings to the standard case of a primordial spectrum with power- law tilt n s , we find a 2a 
detection of strings: fio — 0.11 ± 0.05, where flO is the fractional contribution made by strings in 
the temperature power spectrum (at £ — 10). CMB data give moderate preference to the model 
n s — 1 with cosmic strings over the standard zero-strings model with variable tilt. When additional 
non-CMB data are incorporated, the two models become on a par. With variable n B and these extra 
data, we find that /io < 0.11, which corresponds to G/i < 0.7 x 10 -6 (where fj, is the string tension 
and G is the gravitational constant). 



Introduction. — The inflationary paradigm is success- 
ful in providing a match to measurements of the cosmic 
microwave background (CMB) radiation, and it appears 
that any successful theory of high energy physics must 
be able to incorporate inflation. While ad hoc single-field 
inflation can provide a match to the data, more theoret- 
ically motivated models commonly predict the existence 
of cosmic strings [ljj. These strings are prevalent in su- 
persymmetric D- and F-term hybrid inflation models (see 
eg. 01) and occur frequently in grand- unified theories 
(GUTs) 0. String theory can also yield strings of cos- 
mic extent Hence the observational consequences of 
cosmic strings are important, including their sourcing of 
additional anisotropies in the CMB radiation. 

In this letter we present a multi-parameter fit to CMB 
data for models incorporating cosmic strings. It is the 
first such analysis to use simulations of a fully dynamical 
network of local cosmic strings, and the first to incorpo- 
rate their microphysics with a field theory @, 0| . It yields 
conclusions which differ in significant detail from previ- 
ous analyses based upon simplified models: we find that 
the CMB data Q moderately favor a 10% contribution 
from strings to the temperature power spectrum mea- 
sured at mulitpole £ = 10 with a correponding spectral 
index of primordial scalar perturbations n s ~ 1. There 
are also important implications for models of inflation 
with blue power spectra (n s > 1). These are disfavoured 
by CMB data under the concordance model (power-law 
ACDM which gives n s = 0.951+g ||) and previous 
work seemed to show that this remains largely the case 
even if cosmic strings are allowed (n s = 0.964 ± 0.019 
0). However with our more complete CMB calcula- 
tions, we find that the CMB puts no pressure on such 
models if they produce cosmic strings. Our conclusions 
are slightly modified when additional non-CMB data are 
included, with the preference for strings then reduced. 

CMB calculations. — In the combined inflation plus 
strings case, inflation creates primordial perturbations 



which evolve passively until today but, in the interven- 
ing time period, cosmic strings actively source additional 
perturbations. Given the small size of the observed CMB 
anisotropies, the perturbations may be treated linearly 
and any coupling between those seeded by the two mech- 
anisms can be ignored. The string and inflation per- 
turbations can therefore be evolved via separate calcu- 
lations, yielding two contributions to the CMB power 
spectrum that are statistically independent and so are 
simply added together to give the total power spectrum. 

Calculating the cosmic string component presents a 
challenge because their evolution is non-linear and the 
string width is very much smaller than their separation 
at times of importance for CMB calculations. Previous 
comparisons of the string CMB power spectrum against 
data have relied upon models which neglect the width, 
representing local strings as ID objects and then either 
evolving them according to the Nambu-Goto equations 
appropriate for a relativistic string [10| or employing an 
unconnected segment model (USM) [J ll[ . These USMs 
involve ensembles of unconnected string segments with 
stochastic velocities and with segments removed to mimic 
the time dependence of the string density seen in simu- 
lations. A third approach is to simulate instead global 
strings, which do not localize their energy into the string 
cores. The cores may be left unresolved and field-based 
CMB calculations [T3| have been used elsewhere 13. IH]. 

In we used a field-based approach for local strings, 
via the Abelian Higgs model. We were able to resolve 
the cores and to reach a string separation of ~ 100 times 
their width, which we carefully checked to be sufficient 
to reach a scaling regime. This regime, in which the sta- 
tistical properties of the network scale with the horizon 
size, is of critical importance as it enables the statistical 
results to be applied to the later times required in CMB 
calculations. A great advantage of the field theory is 
that it naturally includes the decay of the string network 
into Higgs and gauge radiation, and the resulting back- 
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reaction on the network. Thus our CMB calculations for 
strings are the first to include a consistent mechanism for 
decay and backreaction. 

A feature of field theory simulations is a very low den- 
sity of string loops [l5j |. in sharp distinction to Nambu- 
Goto simulations on which the conventional cosmic string 
scenario is based. Further work is needed to understand 
the origin of the difference, on which bounds from cosmic 
rays fig or gravitational wave production [j| sensitively 
depend, but CMB calculations depend on the large-scale 
properties, about which there is broad agreement. In- 
deed the USM has enough flexibility to approximate our 
power spectrum: the left hand graph of Fig. 1 of the er- 
ratum to H| is similar to Fig. 13 of [a]. However, the 
USM docs not reproduce the detailed shape of the power 
spectra, nor can it give limits on the string tension fi 
without reference to simulations such as ours. Our cal- 
culations represent a significant step forward in reliability 
and accuracy, deserving careful comparison to the data. 

Data fitting approach. — The form of the cosmic string 
contribution to the temperature power spectrum is shown 
in Fig. [I] where it is compared to observational data and 
the best-fit standard inflation model. The normalization 
of the inflation and string power spectra components are 
free parameters, with that for strings being proportional 
to (Gfi) 2 (where G is the gravitational constant and fi is 
the string tension). For Fig. [T]the normalization of the 
string component has been set to match the data at mul- 
tipole £ = 10, corresponding to G[i = (2.04±0.13) x 1(T 6 , 
a factor of 2-3 higher than the corresponding value from 
previous work [l0, U , [3] ■ Clearly a string component 
this large is ruled out and we hence introduce the param- 
eter /io, the fractional contribution from cosmic strings 
to the temperature power spectrum at £ = 10. 

Recalculating the inflationary component at a particu- 
lar cosmology takes only a few seconds, but for the string 
contribution this takes many hours and it therefore ap- 
pears that a full Markov chain Monte Carlo (MCMC) 
multi-parameter fit is unfeasible. However, following [l7| . 
we fix the form of the string component and vary only 
its normalization, via G\x. Given that any changes in the 
cosmological parameters are small and that the strings 
are sub-dominant, this amounts to a small error in the 
total inflation plus strings prediction, below the uncer- 
tainties in the CMB data [23| and the MCMC results 
are unaffected. We hence use a version of the standard 
COSMOMC [3] code, modified to incorporate the fixed- 
form cosmic string component. 

We primarily consider four different models: two pa- 
rameterizations of the primordial power spectrum, both 
with and without strings. We always allow for variations 
in the Hubble parameter h, the physical baryon and total 
matter densities tl^h 2 and f2 m /i 2 , as well as the optical 
depth to last scattering r. We then either take Harrison- 
Zeldovich (scale-invariant) adiabatic primordial pertur- 
bations with amplitude A s or add the additional freedom 
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FIG. 1: The temperature power spectrum contribution from 
cosmic strings, normalized to match the WMAP data at 
I = 10, as well as the best-fit cases from inflation only (model 
PL) and inflation plus strings (PL+S). These are compared 
to the WMAP and BOOMERANG data. The lower plot is 
a repeat but with the best-fit inflation case subtracted, high- 
lighting the deviations between the predictions and the data. 
Note that the string contribution is identical to that shown in 

Fig. 14 of a , but here has a linear horizontal axis for £ > 100. 



of a power-law tilt n s : A 2 — > A 2 (k/k ) na . This yields 
the two zero-string models which we label as HZ and PL 
respectively, with PL being the established inflationary 
concordance model and HZ being a restriction of this: 
n s = 1. We add strings to these two models yielding 
models HZ+S and PL+S, which therefore have the ex- 
tra parameter (G[i) 2 . Then, in the later stages of our 
discussion, we also consider primordial tensor perturba- 
tions and a finite running of the scalar spectral index 
dn s /d In k, but we will assume negligible neutrino mass 
and flat space throughout. 

Results using only CMB data. — The results when 
using measurements from the WMAP, ACBAR, 
BOOMERANG, CBI and VSA projects Q are illustrated 
in Fig. [5] This shows the marginalized 2D likelihood sur- 
faces for /io versus h, Qbh 2 , A 2 and n s for both HZ+S 
(points) and PL+S (contours). For PL+S, there is a 
significant degeneracy, involving primarily these five pa- 
rameters, that allows large values of /io to fit the data 
The result is f w = 0.11 ± 0.05, which is a 2a de- 
tection of strings. It also yields n s = 1.01 ± 0.04 which 
is significantly larger than in model PL, or the result of 
n s = 0.964 + 0.019 found in 0] for PL+S using the USM. 

FigureQJlower) shows the deviations between the best- 
fit PL+S case, the best-fit PL case and the CMB data. 
Given that the best-fit PL+S case is given by /io = 0.099 
and n s = 1.00 (see endnote [29| for the other parameter 
values), it is clear that not only is n s = 1 under no pres- 
sure if cosmic strings are included, but it is able to fit the 
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FIG. 2: The 2D marginalized likelihood distributions from 
CMB data (only) for /io versus h, Qbh 2 , A 2 and n s . Contours 
show the 68 and 95% confidence regions for model PL+S while 
the 400 MCMC points indicate the prefered region for HZ+S. 
The vertical lines on the h and Q.hh 2 plots show the 68 and 
95% confidence limits from the HKP and BBN measurements. 
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ID 


no. 
par am. 


CMB only 


CMB+HKP+BBN 


Axoff 


evidence 


AXeff 


evidence 


HZ 


5 


+7.7 


0.35 ± 0.03 


+10 


0.120 + 0.009 


PL 


6 





1 





1 


HZ+S 


6 


-3.9 


7.3 + 1.2 


+0.9 


0.68 + 0.12 


PL+S 


7 


-3.9 


1.2 + 0.1 


-1.6 


0.19 + 0.01 



TABLE I: The Axoff an d relative Bayesian evidence values for 
the examined models using the CMB, HKP and BBN data. 



data moderately better than the n s = 0.952 best-fit under 
model PL. Indeed, when the maximum likelihood values 
£ max are compared via Axl s = -2 ln(£^+ s /£^ x ), we 
obtain A%^ ff = —3.9 at the expense of a single extra 
parameter. However, as the PL+S best-fit value of n s 
is extremely close to one, HZ+S has an almost identical 
£ max value. Therefore model HZ+S gives Ax^ = —3.9 
relative to the concordance model with zero cost in terms 
of the number of parameters. 

A more complete analysis of the freedom in a model is 
provided by its Bayesian evidence value [3] and Liddle 
et al. [h| have previously used this statistic to demon- 
strate that WMAP data does not actually rule out model 
HZ, despite the n s = 0.951^;^ result returned under 
the standard model PL. Here, we calculate evidence ra- 
tios for our four models using the Savage-Dickey method 
[H with flat priors of < / l0 < 1 and 0.75 < n s < 1.25, 
giving the results shown in the table. We find that the 
relative evidence of PL+S to PL is barely distinguish- 



able from unity, as expected for merely a 2a detection 
of strings. However, model HZ+S has a Baycs factor 
of 7.3 ±1.2 relative to PL and is therefore moderately 
preferred. That is, a finite string component is favored 
by CMB data over a tilted power spectrum and the re- 
sult: /io = 0.10 ± 0.03 from model HZ+S is therefore of 
interest. 

Use of non-CMB data. — We must also check that 
these conclusions remain valid when non-CMB data are 
included and we hence consider that the Hubble Key 
Project (HKP) yielded h = 0.72+0.08 (H)]. Further, mea- 
surements of deuterium abundance in high redshift gas 
clouds, combined with big bang nucleosynthesis (BBN) 
calculations, gives n b h 2 = 0.0214 + 0.0020 [H[ and while 
similar determinations using other light isotopes do not 
yield global concordance, it is still interesting to consider 
this measurement also. Figure [2] shows these two mea- 
surements via vertical lines in the relevant plots and it is 
clear that they each disfavor large values of /io in model 
PL+S. It is also evident that they lower the preference 
for model HZ+S since the majority of the plotted MCMC 
plots lie at least la from these two results. 

With these data included, model PL+S now yields 
/io = 0.05±g g| or /io < 0.11 (95% confidence) and 
the 2a detection is removed. However, the result of 
n s = 0.97 ± 0.02 still does not ride out n s > 1 with 
any confidence (cf. n s = 0.953 ±0.015 obtained using the 
USM with these data 0). The Bayes factor for model 
HZ+S relative to PL is reduced, but only to 0.68 ± 0.12, 
leaving HZ+S on par with the standard model. 

We also incorporate galaxy survey data via the matter 
power spectrum, although there are uncertainties over 
the use of such data when strings (or other defects) are 
included [l7j]. However the CMB constraints, together 
with our calculated string contribution to the matter 
spectrum, imply that strings make a negligible contri- 
bution to the matter power spectrum on large scales (as 
is also the case using USM calculations ll|). These are 
the same scales where the zero-string case needs no cor- 
rections for non-linearity, which have been questioned in 
[23I ]. We therefore conservatively include SDSS Lumi- 
nous Red Galaxy data 0] for only k/h < 0.08 Mpc -1 
finding that it leaves our results essentially unchanged: 
f w = 0.10 ± 0.04 and n s = 1.00 ± 0.03 for PL+S with 
an evidence value of 7.7 ± 0.7 for HZ+S. Including also 
data for 0.08 < k/h < 0.2 Mpc -1 and non-linear correc- 
tions gives /10 < 0.11, n s = 0.97 + 0.02 and evidence 
0.50 ± 0.05 but the use of the non-linear regime makes 
these results less reliable. 

Hence, while we await further updates from the ob- 
servational community regarding these additional data, 
even with them included, model HZ+S remains compet- 
itive relative to PL. 

Tensors and running. — When the freedom for a non- 
zero primordial tensor contribution is incorporated as a 
generalization of model PL, tensor modes give a negligi- 
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ble (and possibly zero) improvement in the fit to CMB 
data. However they do raise the allowed n s to 0.98 ±0.03 
[U (CMB only), which is a greater effect than the USM 
strings of Q. As an addition to PL+S, tensors are more 
preferred but again they increase the allowed n s values. 
For the CMB+HKP+BBN case we find n s = 0.99 ±0.02, 
hence even the BBN data puts no pressure at all on 
n s > 1 when both strings and tensors are included. 

Adding finite dn s /d\nk (running) to model PL does 
give a marginal improvement to the fit, with CMB data 
preferring a slight negative running Q . This lowers small 
and large scales relative to intermediate ones and may 
hence be thought to have a similar effect as strings. How- 
ever, adding strings smooths out the acoustic peaks and 
in fact there is little correlation between dn s /d\ak and 
/io- Hence we find that the above results are barely af- 
fected by finite running. 

Conclusion. — By including cosmic strings, we find a 6 
parameter model with n s = 1 that performs better than, 
or about as well as, the established concordance model, 
and that the latest data does not necessarily favor n s < 1. 
We also find that, when incorporating the (debatable) 
deuterium BBN result, the cosmic string contribution is 
constrained to /io < 0.11 or Gfi < 0.7 x 10~ 6 . Even 
at this level it is likely that cosmic strings will be soon 
detectable using the B-mode polarization of the CMB Q 
and we await future data releases with great excitement. 

Finally, we note that our bounds have been derived 
only for classical Abelian Higgs strings with equal vector 
and scalar particle masses [5j, and that, for example, F- 
term inflation may be more accurately treated using sim- 
ulations with different values. Similarly, different CMB 
predictions for strings may be found with (p, g)-string 
networks from string theory or from other models 
such as semi-local strings |26j . A confirmed string detec- 
tion would open up the challenge of differentiating the 
models and hence learning a great deal about inflation 
and high energy physics. 
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The absolute uncertainties in the string power spectrum 
calculations are also smaller than those of the data. 
Recalculating the string component at the parameters re- 
quired for high and low /io yields an insignificant change, 
confirming the validity of the method employed. 
The remaining best-fit values are h = 0.82, Q^h 2 = 
0.0255, n m h 2 = 0.123, r = 0.11, and A 2 = 20 x 10~ 10 . 



